A necklace of Wulff shapes 
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Abstract: In a probabilistic model of a film over a disordered substrate, Monte-Carlo simulations show 
that the film hangs from peaks of the substrate. The film profile is well approximated by a necklace of 
Wulff shapes. Such a necklace can be obtained as the infimum of a collection of Wulff shapes resting on 
the substrate. When the random substrate is given by iid heights with exponential distribution, we prove 
estimates on the probability density of the resulting peaks, at small density. 
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1. Introduction 

A problem in the science of coating is to characterize the surface of a coated material as 
function of the substrate surface and properties of the coating material (polymer, resin, 
metal. . . ). The topography of a substrate has an important influence on the properties of 
the considered material in terms of lubrication, optical properties, wetting... Moreover, it 
is often desirable to coat this substrate with a thin film to protect the material. Of course, 
the topography of the thin film surface and of the substrate do not have to be the same. 
This will be a function of the film thickness. 

A statistical mechanical model of a film requires at least two parameters, one asso- 
ciated with a surface tension, the other with a pressure difference or chemical potential 
controlling the film thickness. The version of the Solid-On-Solid model introduced by 
Abraham and Smith [AS1, AS2] is the simplest such model. In Section 2, we use it first 
to model a disordered substrate and then, with a different set of parameters, to model a 
film on top of the substrate. Numerical simulations show that typical configurations look 
like a necklace of Wulff shapes suspended from the peaks of the substrate. This motivates 
Section 3, where only the substrate is random, and the film is defined as the infimum of 
a collection of Wulff shapes over the substrate. In Section 4 we prove estimates on the 
density of relevant substrate peaks, when the pressure difference goes to zero. In section 
5 we give a Gibbs formulation for the probability density of substrate peak localizations 
and heights. 

2. Solid-On-Solid film over quenched Solid-On-Solid substrate 

The substrate is a one-dimensional Solid-On-Solid model with Hamitonian 

H X = J X \h\-h)\ + K x Y,h\ (2.1) 

l«-il=i 
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where h\ e R+ is the height of the substrate surface at point i e Z. It can also represent 
the height of a first coating, fixed before the later deposit of a film. 

The substrate being generated according to the Gibbs measure with Hamiltonian (2.1) 
at temperature kT = 1, and then quenched, a film is deposited and thermalized according 
to the Gibbs measure with Hamiltonian 



Ho 



J 2 \h 2 i-h 2 j\ + K 2 J2hi 



(2.2) 



I*-j'I=i 



at temperature kT = 1. The film height hf G IR+ at point i obeys hf > h\. This is a 
grand canonical ensemble where the film volume is controlled by K 2 while being allowed 
to fluctuate. In the thermodynamic limit, the properties will be the same as if obtained 
from a fixed volume ensemble where ^{h 2 — h\) is fixed. 

The resulting model of a film on a quenched substrate is studied by Monte-Carlo 
simulation with periodic boundary conditions and a heat bath algorithm. Fig. 1 shows 
one substrate h 1 and, on top, thermal averaged films h 2 at various values of the pressure 
K 2 . 
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Fig. 1: Substrate J\ = 1, K\ = 0.5; film J 2 = 30 and from top to bottom K 2 = 0.2, 0.5, 1, 2. 

For small enough K 2 /J 2: the film appears to hang as from a set of telegraph poles of 
random heights. Between two successive poles, the film profile can be checked to be very 
near a Wulff shape associated with (2.2) at the corresponding value of K 2 . Indeed Fig. 2 
shows a portion of substrate, a thermal averaged film with J 2 = 5 and K 2 = 0.125, and a 
necklace of Wulff shapes: Each piece is a translate of one and same Wulff shape, with the 
same J 2 and K 2 as the film. The parametric equations of the Wulff shape are [BN,DD]: 

x(tan^) = — -^r ^ ^(tanfl) 

. (2-3) 
-cr(tan0)J 



K 2 d(tan0) 
r(tan#) = — — - (a{ia,n9) — tan 9 



d 



d(tan0) 
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where a(tan6>) is the projected surface tension, which for the Solid-On-Solid model takes 
the form [DD] 

<r(tan0) = /(tan0) - log ^^ + 2 ) (2.4) 

with 

/(tan0) = (l + (J 2 tim6) 2 ) 1/2 -1 (2.5) 
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Fig. 2: Substrate J x = 0, K x = 0.5; film (solid line) and Wulff shapes J 2 = 5, K 2 = 0.125. 

The Wulff necklace is clearly a good approximation, except when the substrate rises 
near but below the film surface, like near x ~ 80 or x ~ 140 on Fig. 2, which causes some 
entropic repulsion. As K 2 decreases, the number of substrate peaks visible from the look 
of the film decreases (see Fig. 1). The relevant peaks are large deviations events of the 
substrate disorder. But there is no obvious a priori criterion to tell which h\ will emerge, 
depending upon K 2 . The film tension leads to correlations between the poles. 

In the following sections, we simplify the model in order to understand the selection by 
the film of substrate peaks as K 2 varies. When the substrate has short range correlations, 
its correlation length may be taken as a basic unit of length. A natural simplification is to 
assume that it is also the lattice unit, and that the substrate is i.i.d. on this scale: J\ = 0. 
Then, by scaling the unit of height, there is no loss of generality in taking K\ = 1, leaving 
two independent parameters J 2 and K 2 . 

The second and main simplification is to consider that when a Wulff shape between 
two peaks hangs strictly above the substrate between the two peaks, then it is not affected 
by the substrate between the two peaks. Some fluctuation effects (entropic repulsion) are 
neglected here, an approximation which is better justified when J 2 is large. 

Our aim will be to estimate the density of relevant peaks as function of K 2 , with 
different kinds of Wulff shapes. And further to get an idea of the probability distribution 
of peak localizations and heights. 
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3. A necklace of Wulff shapes 

Let hi, i G Z be an iid sequence of exponentially distributed random variables of mean 1, 
and let 

S = {(i, z) G Z x R : z < hi) (3.1) 

be the random set representing the substrate. The upper index 1 for the substrate height 
h\ = hi is now omitted. 

Let 1 < a < oo and let W : ] — a, a[— > R be a continuous even function, strictly 
increasing on [0, a[, with W(0) = 0. If a = oo, we require W(x) > — b + X\x\ a for some 
b, A, a > 0, for all x. Examples: 



Cone : 


W(x) 


= X\x\ 


Parabola : 


W(x) 


= Xx 2 


Semi — circle : 


W(x) 


= x~' 


SOS Wulff shape (z(x) = (2.3)) : 


W{x) 


= z(x) 



(3.2) 



1 

K~2 



2 

T 2 



with a = X 1 for the circle and a = J 2 /K 2 for the SOS Wulff shape. Translating the graph 
of such a function by (x*, h*) G R 2 will define W(x*, h*; •), so that W(0, 0; •) = W(-) and 



W(x*, h*; x) = h* + W(x — x*) 



(3.3) 



The graph of this function, as a subset of R 2 , is denoted W{x ie ,h if ). The film over the 
substrate (or the coating of the substrate) is then defined almost surely as the graph / of 
the function 7(x) defined in the following proposition: 

Proposition 3.1: Under the hypotheses on W(-) stated before (3.2), let 



I(x) = inf jw^x*, x) : (x*, h*) G R 2 , |x* — x\ < a, W(x*, /i*) fl 5 = j 



(3.4) 



whenever the inhmum is taken over a non-empty family, and I(x) = +oo ifW(x*, h*)(lS ^ 
V(x#, h*). Then, almost surely, the inhmum is attained and I(x) < oo Vx G R. 

Proof. Straightforward. The hypothesis on W(-) could be weakened to W(x) > —b + 
Alog|x| for suitable A. In the following, we shall instead strengthen the hypothesis to 
make W convex. 
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Fig. 3: 7(x) with W(x) = A|x|. 
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Fig. 4: /(x) with W(x) semi-circular. 

In words: We have a "Wulff shape", symmetric about a vertical axis. Above each 
x* G M, a "Wulff shape" is taken down from +oo, until it touches the substrate. The 
film is the infimum of the resulting collection of Wulff shapes. The film height I(j) at an 
integer point j models the thermal average of h 1 in the preceding section. 

A Wulff shape in one dimension is a solution to the second order differential equation 
for the function W(x), 

4- 7 d , r a(dW/dx) =K 2 (3.5) 
dxd{dW/dx) K 1 ' y ' 

where cx(-) is the projected surface tension, or interface free energy per unit length of 
interface projected onto the x-axis, as function of the slope dW/dx. The parameter K 2l 
conjugate to the film volume ^2hf in (2.2), is the pressure, or pressure difference with 
the medium above, Ap. The surface tension should have convexity properties such that 
the solution W(x) to (3.5) is a convex function, in fact typically strictly convex, which we 
assume henceforth, except for the special case W(x) = X\x\. 

A Hamiltonian of the form (2.2) but quadratic in (hj — h%) gives 

a [dW/ dx) = J 2 [dW/ dx) + constant 

so that a solution to (3.5) in this case is W(x) = Xx 2 with A = K^j^Ji). Solutions to the 
Wulff equation (3.5) generally scale as 

W K2 {x)=K- 1 W 1 {K 2 x) (3.6) 

where W\{x) is a solution to (3.5) with K 2 = 1. The semi-circular shape in (3.2) follows 
this scaling, with A proportional to K 2 . 

When a Wulff shape is taken down from infinity above x = until it touches the 
substrate *S, say at (jo, hj ) with j Q < for definiteness, it stops as 

w(o, MO); x) = h jo - W(j ) + W(x) (3.7) 

We can then let it slide to the right by an amount x*, keeping contact with (jo, hj ), as 

W(x*, h*(x*);x) = h jo - W(j - x*) + W(x - x*) (3.8) 
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At a given x > 0, this is a strictly decreasing function of so long as x* < x. It is also a 
strictly decreasing function of x* for x* > x, provided W(-) is a strictly convex function: 
For > x*, 

W(y*) - W(y* - x) > W(x m )-W(x m -x) 

X X 

So we let the Wulff shape slide to the right keeping contact with (jo, hj ) until it touches 
a second point (ko, hk ): 

hk - h jo = W(k - x,) - W(j - x,) (3.10) 

We thus get the unique Wulff shape going through (j , hj ) and (ko, hk ), which we denote 
W(jo, hj , ko, hk ; x). If the first contact point is (ko, hk ) with ko > 0, instead of (jo, hj ) 
with jo < 0, the same construction works symmetrically, sliding to the left. In any case 
we get a Wulff shape at a minimal height, so that, almost surely, 

I(x) = W(jo,h jo ,ko,h ko ;x), jo < x < k (3.11) 

The proposition below follows easily: 

Proposition 3.2: Let W : R — > R be a continuous even function, strictly convex with 
W(0) = 0, or W(x) = X\x\ with A > 0. Then I(x) defined in (3.4)(3.1) is also, almost 
surely, 

I(x) = supjwC?, hj, k, h k ; x) : j, k e Z, j < x < /c| (3.12) 

Proof. Starting from (3.4), we have I(x) = W(x*, h*; x). This must be also I(x) = 
W(j,hj,k,hk',x) for some j, k: Otherwise (3.7-11) would give a smaller I(x). And this 
W(j, hj, k, hk] x) must be the supremum, otherwise it would intersect S (as shown in more 
detail in the proof of Lemma 5.1 below). 

Remark: For the first two examples in (3.2), we have respectively 

W(j, hj, k, hk; x) = max j/ij — X(x — j), hk + X(x — fc) j (3.13a) 

W(j, hj, k, h k ;x) = X(x - j)(x - k) + 7 — - hj + ] - — % h k (3.136) 

k — 3 j - k 

The random interface I(x) also defines, as a marginal, a point process of interest: 
Proposition 3.3: Under the same hypotheses as in Proposition 3.2, let 

B = /n<S= {ieZ:/(i) = AiJ (3.14) 

Then 

B = {i EZ: hi>W(j,hj,k,h k ;i) V j < i < k} (3.15) 

and, almost surely, B can be written as B = {b n } ne x with b n+ i — b n > 1 Vn, and bo = 
min{6 n : b n > 0}. 



Proof. If i belongs to (3.14) then Proposition 3.2 implies that it belongs also to (3.15). If 
i belongs to (3.15), then we start from 

hi = W(i, h { ; i) > W(j, hj,k, h k ;i) V j <i<k 

and slide W(i,hi;i) to the right following (3.7-11) until it touches S at (k,hk) with hi = 
W(i, hi, k, hk', i). Then Proposition 3.2 implies that i belongs also to (3.14). 

The set B is also the set of points which can be obtained like j or k in (3.5)-(3.8), 
starting from any x e R, not just x = 0. 

4. Estimates 

Proposition 4.1: Let W{x) = X\x\. Then 

Urn = i (4.1) 

A\o A 2 v 1 

Proof. 

P(0 e B) = P(\) > max{hj - X(x - j), h k + X(x - k)} V j < < fc) 
= F(ho>hi-X\i\ VzeZ*) 

roc oo 2 

= / dxe- x Y[(l-e- x - Xl ) 

J ° i=i K 1 

— Upper bound: 



(4.2) 



oo 2 



i=l 

so that 



- e-*- Xi ) < e~ 2 ^ Sr=x = e" 2e "^ (4.3) 



/" 1e~' e ~ A l-e" A / 2 e ' A \ 

P(0GB)< rfxe- x e" 2e = l—h-e" 2 !^) (4.4) 

Jo 2e V / 

Lower bound: For x > 1, we can use 1 — e > e _<E_e2 , with e = e _:E_A \ Then 

1 - e" x - A M >e" 2e ^=i e " 2e ^^i e =e ^t^x 2e T -^ x ^ 



so that 



/oo 
<fcce _ V 



2e~ x e \ -2e' 



/AiniA 1 -A „ -2A 

^ e -^Tt^x-^Tt^x (4.6) 



2c 
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The upper and lower bounds together imply (4.1). 

Proposition 4.2: Let W{x) = Xx 2 . Then 3 a, b > such that VO < A < 1/4 

a R ^ Ty <P(0GB)<3^ + 6A (4.7) 

Remark: In view of the proof of the upper bound, see below, we expect that P(0 G B) is 
of order \/X for small A. 

Proof. 

F{OeM) = F(h >W{-j,h_ j ;k,h k ;0) V j, k > l) (4.8) 

with 

W(-j, h-j-k, h k - 0) = — - h-j + —J— h k - Xjk (4.9) 

— Upper bound: 



(4.10) 



P(0 £ M\h = x) = P( 3 j, k> : x < W(-j, h_f, k, h k ; 0) ) 

> P( 3 j, k > : h-j > x + Xj 2 , h k > x + Xk 2 ) 

oo 2 

= {l-l[¥(h. J >x + Xj 2 )) 

3 = 1 

where we used 

W(-j,x + Xj 2 ;k,x + Xk 2 ;0) > x (4.11) 
and W(- ■ •) increasing in h-j and in hk- Then 

l[F(h- J >x + Xj 2 ) = l[(l-e- x - x i 2 )<e- e ~~ X ^ e ~ * < e^^-D (4.12) 
j=i i=i 

so that 

P(0 6 B) < J°° dxe-'(2e-'~"<-iVJ-» - e -^"(ivT-i>) 

= J r i *( 2e -.(»v¥-»- e ->.(*^)) (413) 



Lower bound: From the Harris-FKG inequality [H,FKG], 

oo 

P(Vj,fc>0: x>W{-j,h-y,kM^)) > Y[F(\/k>0: x > W(—j, h-j-,k, hk] 0) ) 

i=i 

(4.14) 



poo ^~ 

P( VA; > : x > W(-j, h-y, k, h k ; 0) ) = / dy e~ y TT P( x > W(-j, y; k, h k ; 0) ) 

Jo k=i 

r x+Xj 2 oo 

/ dye-y\[(l-e- Xk2 - k -^ 2 -y)-^ 
^° k=i 

„ x -\-\j 2 oo oo 

^° fc=i fc=i 



> 



_ _ 2 _ 

where we used 1 — e > e e 6 , with e < e x and x > 1 henceforth. Then 



(4.15) 



oo oo 



-X 



fc=i fc=i e — 3 1 ^ 

F(Vk>Q:x>W(-j,h- j ;k,h k ;Q))> * dye~ y ~ »+**-v 

Jo 

~7I 



g x + \j z / _ x + je 

> — 



1 + 



2je- 



2 



t+ZF^j (4.17) 



/ je X lie X J! + Aj 2 je x ( , y -2n 2je 31 \ 

>exp(-^— J A . 2x2 -e T (x+Aj } ^) 

V x + Aj 2 (x + Aj 2 ) 2 / 

where (1 — ay) -1 < 1 + 2ay for < ay < 1/2 has been used before integration over y, and 
then again 1 — e > e~ e_e and also 1/(1 + X) > e~ x . Now 



A" 1 



x + A 1 e x 



|TP(VA;>0: x>^(-i,/i_ j; /c,/i fc ;0)) >exp(-^— In , 

+ VAX (4.18) 

-e-/ 2 - - 2 




2A(x + A) V A V A 



The sums over j e Z + in the exponential were bounded by corresponding integrals over 
[1,A _1 ] or R_|_ plus a bound of the maximum of the integrand. Then with A < 1/4 and 

x > 4, 



A' 1 



J]P(Vfc>0: x>Ty(-j,/i- j; fc,/i fc ;0)) >exp(-— lnA- 1 --^- V - 2 j 

(4.19) 

For j > A -1 , the range of integration in (4.15) is chosen as < y < x + Xj 2 — j -In 2. 
Then only the first few k = 1, 2, . . . play a role for the event 3 k : x < W(—j, y; k, h k ', 0). 
Precisely: 

f; e -^-k(x + Xf-y)- X < __g < e -x-^±M^ + 2e -,-2^M^ (4 _ 2Q) 

fc=i e j — 1 



x+\j 2 -j In 2 



( x x + x i 2 -y _ 2 X+X: > 2 ~ v \ 

ay exp I — y — e — 2e i ) 



,x+Aj 2 -jln2 , 



1 ' ' (4.21- 

_ x _£+m£. 

= 1 _ e -(x+Aj 2 -jln2) 3 e (* c -(x+A. 7 2 -, 7 ln2)(l-l/, 7 )j 

> 1 - 5e" x - Aj > exp(-5e- x - Aj - 25e~ 2x - 2Xj ) 
for A < 1/4 and x > 2. Then 

f j P( VA; > : x > W(-j, h-j\ k, h k ; 0) ) > jj exp(-5 e~ x ~ X3 - 25 e - 2x ~ 2Xj ) 

j>A-i j>A-i 

Se" 1 " 1 25e" 2x - 2 \ 



/ be 



e -A ^ _ e -2A 



Putting together (4.19) and (4.22) and integrating with d(e x ) over an interval 



(4.22) 



yields the lower bound and completes the proof of the proposition. 
5. Gibbs measure 

So far the equivalent definitions of B, through (3.14) or (3.15), require a knowledge of the 
whole system in order to decide whether a point i G B. Here we will define local Gibbs 
factors where (b n ,hb n ) is coupled to (b n -i, ht >n _ 1 ) and (b n +i, ht, n+1 ) only. We first derive 
an equivalent definition of B: 

Proposition 5.1: B defined in Proposition (3.3) is almost surely the only countable 
ordered subset ofZ, denoted B = {& n } n eZ> with bo = min{6 n > 0}, obeying the following 
two conditions: 

hi <W(b n ,h bn ,b n+1 ,h bn+1 ;i) V b n < i < b n+1 (5.1) 

h n > W(b n -i,h bn _ 1 ,b n+1 ,h bn+1 ;b n ) V neZ (5.2) 

B is the minimal subset of Z such that the collection of Wulff shapes suspended from this 
subset lies above the rest of the substrate, i.e. the minimal subset ofZ satisfying (5.1). 

Proof. From Propositions (3.2) and (3.3), B is almost surely the only countable ordered 
subset of Z satisfying (5.1) and 

h bn >W(b n - p ,h bn _ p ,b n+q ,h bn+q ;b n ), V neZ,p>l,q>l (5.2') 

We only need to prove that (5.2) implies (5.2'). Let us assume two instances of (5.2), 

h bl >W(h bo ,h b2 ;h) (5.3a) 
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h b2 >W{h bl ,h bs ;b 2 ) (5.36) 
and, contradicting (5.2') with p = 2, q = 1, 

h b2 <W(h bo ,h b3 ;b 2 ) ??? (5.3') 

Let us show that (5.3a — 6) with the absurd (5.3') would imply that W(h bl , h b3 ; •) and 
W(hb , hb 2 ] •) have two intersections. Indeed 

W(h bl , h ba ; 6i) > W(h bo , a &2 ; h) using (5.3a) 
VT(/i fcl , h ba ; 6 2 ) < W^(/i 6o , h b2 ; 62) using (5.36) 
W(h bl , a b3 ; 63) = W(h bo , a b3 ; 63) 

= W(h bo ,W(hb ,hb 3 ;b 2 );bs) 
>W(h b() ,h b2 ;b 3 ) using (5.3') 

so that W(h bl , hb 3 ; •) is higher than W(h bo , h b2 ; •) at 61 and 63 and lower at 62, which 
implies two intersections, impossible for two Wulff shapes which are translates of one 
another. Therefore (5.3') cannot hold true; (5.2) with p = 1, q = 1 implies (5.2) also with 
p = 2, q = 1. The argument extends easily to all p, q and the proof of Proposition (5.1) is 
readily completed. 

Let l n = b n — 6 n _i and x n = hb n . A Gibbs measure formulation for {I„,i n } ne 2 starts 
from i.i.d. a priori measures: Counting measure on Z + for each Z n , exponential distribution 
exp(— x n )dx n on R + for each x n . And a product of local Gibbs factors, 

Y[F(x n J n+ i,x n+ i) Y\_G(x n -i,l n ,x n ,l n+ i,x n+ i) (5.4) 

n n 

with 

h-i 

F(x , luxt) = Y[ (l - e -w(o,* ,li,*i;i)^ (5.5) 
i=i 

and 

G(x_i, Z , X , /l, Xi) = l Xo > VK(-/ o ,a;-i,«i,a;i;0) (5-6) 

Let us define explicitly a finite volume Gibbs measure with free boundary conditions: 

Proposition 5.2: Let hi, i G Z be an iid sequence of exponentially distributed random 
variables of mean 1. Let W : R — > R be a continuous even function, strictly convex with 
W(0) = 0. Let 

I[0,L](i) = sup{^(j, hj, k, h k ; i) : < j < i < k < L} (5.7) 

and 

B [0) l] = |i G [0, L] n Z : 7 [0iL] (i) = = {6 , . . . , b N } (5.8) 
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with bo = 0, &at = L, and N > 1. Let Z n = 6 n — 6 n _i and x n = ht, n - Then 
f^N, {/!,..., /^r}, {xo, . . . , is distributed according to 

N N-1 N-1 

) (5-9) 

1 

where the partition function is 

L r N N-1 N-1 

S [0,L] = / II rf;rne ~ X " II F ( x nJn+l,X n+1 ) JJ G{x n _ 1 , l n , X n , l n+1 , X n+1 ) 

N=lh...l N ^ 1 

(5.10) 

and the sum over the positive integers l\ . . . In is constrained by l\ + . . . + In = L. 

Proof. Proposition 5.2 is a simple corollary of Proposition 5.1. 

Such a model is solvable in principle. A natural first step is to change to a pressure 
ensemble with exp(— pL) and L random in order to get rid of the global constraint over 
h . . . In- 
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